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, ^, Abstract. We present a new model of the diffuse Galactic synchrotron radiation, concentrat- 

ing on its angular anisotropics. While previous studies have focussed on either the variation 
^ of the emissivity on large (kpc) scales, or on fluctuations due to MHD turbulence in the 

QQ interstellar medium, we unify these approaches to match the angular power spectrum. We 

t — note that the usual turbulence cascade calculation ignores spatial correlations at the injection 

^-^ scale due to compression of the interstellar medium by old supernova remnants — the 'radio 

loops' only four of which are visible in radio maps. This new component naturally provides 
Q the otherwise missing power on intermediate and small scales in the all-sky map at 408 MHz. 

Cr^ Our model can enable more reliable subtraction of the synchrotron foreground for studies of 

^ CMB anisotropics or searches for dark matter annihilation. We conclude with some remarks 

^ on the relevance to modelling of the polarised foreground. 

X 
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1 Introduction 

Radio maps of the sky contain a wealth of information on non-thermal astrophysical pro- 
cesses. The sources of this radiation range from compact objects like pulsars with magne- 
tospheres extending over hundreds of kilometers, to radio galaxies hundreds of kiloparsecs 
across, thus spanning over 16 orders of magnitude in size. The total brightness of the radio 
sky is dominated by diffuse emission from the Galaxy at frequencies between tens of MHz 
and a few GHz [1]. It is produced by synchrotron radiation of cosmic ray (CR) electrons and 
positrons^ gyrating in the Galactic magnetic field (GMF) and for typical field strengths of 
(9(/iG), the CR electron spectrum is probed between hundreds of MeV and tens of GeV. The 
frequency and angular dependence of this diffuse emission encodes both the properties of the 
CR electron sources, e.g. their spectrum and spatial distribution, and of the GMF and its 
effect on the propagation of CRs. Therefore, studies of the diffuse radio background can in 



^In the following, when using the term 'electrons' we refer to both electrons and positrons. 



principle determine crucial parameters of the propagation of CRs, e.g. the size of the (radio) 
halo in which they are confined. 

The diffuse synchrotron emission at GHz frequencies also constitutes an important fore- 
ground for studies of cosmologically important angular anisotropics in the cosmic microwave 
background (CMB) [2]. Only in a narrow window around 60 GHz do the CMB temperature 
anisotropics dominate over synchrotron, dust and thermal bremsstrahlung emission from the 
Milky Way. In particular, detection of the elusive B-mode polarisation, a diagnostic of gravi- 
tational waves expected to have been generated during primordial inflation, requires detailed 
understanding and rehable subtraction of the Galactic 'foreground' [3, 4]. This is also a 
limiting factor for the sensitivity of radio searches for dark matter annihilation signals as has 
been illustrated in the context of the 'WMAP haze' [5]. 

As regards the GMF, it is customary to distinguish between an ordered and a 'random' 
component; the former is coherent on (9(kpc) scales while the latter varies on scales up to 
(9(100) pc. Consequently, the specific synchrotron emissivity, i.e. the power emitted per 
unit volume and frequency, varies on two (disparate) scales. On small (< 100 pc) scales 
the variation is almost entirely due to the turbulent nature of the interstellar magnetic 
field. It is believed that from the 'outer scale' L ^ 100 pc on which energy is injected 
into the interstellar medium (ISM) — most likely from supernova remnant (SNR) shock 
waves — the turbulence cascades down to smaller and smaller scales. Both theoretical 
arguments [6, 7] and MHD simulations suggest a power law for the turbulent energy in 
magnetic fields, dB(k)^ /d^k oc k~^, with a = 11/3. The spatial correlations of the magnetic 
field lead to angular correlations of the synchrotron flux J, C('0) = (J(n)J(n')), where n 
and n' are two directions on the sky and ^ = arccos (n • n') is the angle between them. 
A convenient representation for studying these correlations is the angular power spectrum 
(APS): 

Ci = /d(cos^)PKcos^)C(^) , (1.1) 

i.e. the Legendre transform of the angular two-point correlation function C{t/j) — structure 
on angular scales '0 is encoded in multipoles / ^ 7v/t(j. The APS of a sky map J(0, (/)) can 
also be calculated from its spherical harmonics coefficients aim'. 

^^ = ^[^ E 1"^-!' ^^^^^ ai^ = JdnYi*^ie,<i>)Jie,<p)- (1.2) 

m=—l 

This is equivalent to the definition in eq. 1.1 for a (statistically) isotropic sky; here, Yir^iO^ (j)) 
denotes the spherical harmonics function. 

For a power law spectrum of magnetic fluctuations, quasi-linear theory for wave-particle 
interactions predicts a rigidity-dependent diffusion coefficient, which is also a power law, 
Dxx ^ T^^ with 5 = 4 — ^ = 1/3, which seems to be borne out by observations of nuclear 
secondary-to-primary ratios [8] (although [9] reports a different result). Therefore, for the 
energies relevant here, i.e. E < 100 GeV, any structure in the distribution of CR electron 
sources is washed out by diffusion, such that the CR electron distribution is quite smooth on 
scales of < 100 pc. 

On large (> 1 kpc) scales, however, there are variations in both the CR electron dis- 
tribution and the GMF. The variation of the electron density is determined by the spatial 
distribution of CR electron sources and the details of the propagation model, e.g. shape and 
extent of the propagation volume, diffusion coefficients and speed(s) of convective wind(s) 



if any. For the GMF, the ordered component fohows the the Galactic spiral arms, while 
the RMS value of its turbulent component is expected to decrease with distance from the 
Galactic centre and with distance from the disk. (See [10] for a review on (Galactic) magnetic 
fields.) 

Attempts at modelling the diffuse synchrotron background have focussed on the dis- 
tribution of emissivity on either very small or very large scales. It is known [11-14] that 
(assuming a statistically isotropic sky) the power law of the ISM turbulence is reflected in a 
broken power-law for the APS: below some critical multipole /cr? Ci oc /~^, whereas for larger 
multipoles (smaller angular scales): Ci oc l~^. The critical multipole depends on the 'outer 
scale' of turbulence L and the scale height of the turbulent medium, R: /cr ^ {2tiR/L). How- 
ever, comparison with observations is hampered by the fact that the observed synchrotron 
sky is not statistically isotropic, due to large-scale variations of emissivity. Studies of MHD 
turbulence in the APS have therefore constrained themselves to patches of the sky, e.g. by 
cutting out bands around the Galactic plane. Clearly, this cannot describe the APS down 
to small /; in fact, attempts to fit the APS by a broken power law are able to reproduce the 
observed / < 20 behaviour, /(/ + 1)C/ ^ constant, only by assuming an unrealistically small 
scale height: R < Ikpc [15]. 

On large scales, two approaches have been adopted in the literature: Motivated by 
the need for a realistic model of the polarised emission, one set of studies [16-19] employ 
sophisticated models (of the ordered component) of the GMF constrained by detailed analyses 
of rotation measures (RM) and starlight polarisation, but make rather simplistic assumptions 
for the spatial distribution of relativistic electrons. The other set of studies [20-23] start with 
a detailed model for CR propagation, reproducing locally measured spectra of protons, nuclei 
and electrons, and fix the (random component) of the GMF by fitting to the strength of the 
synchrotron radiation. Most of the emission on large scales would appear to be produced by 
the random component — in fact the magnitude of the ordered component as obtained from 
RM data is too weak to explain the radio flux, as has been recognised for some years [24]. 
However, the random component has not been modelled as a true random field but only in 
terms of its RMS value. While this modelling of the large-scale component can reproduce the 
global behaviour, there is mismatch on intermediate scales of tens of degrees. To downplay 
this, angular profiles of the radio background are usually averaged over substantial parts 
of the sky before comparison with data. The remaining discrepancies on smaller scales are 
usually attributed to "local structures" but there have been no attempts to model these. 

It is thus necessary to model the emission consistently on both large and small angular 
scales. There must also be additional contributions at intermediate angular scales from the 
thousands of old supernova remnants (SNRs) which have expanded to large sizes of (9(100) pc 
in the hot tenuous ISM before coohng and turning radiative [25]. Although these old SNRs 
may not actively accelerate CRs, they would compress the GMF in their 'supershells' by 
large factors, and with it the tied CR electrons, which would moreover be then betatron 
accelerated to higher energies [26]. As the synchrotron emissivity scales as the square of 
the magnetic field strength and of the correlated increase in CR electron energy [27], the 
radio emission from such old SNR shells can be significantly amplified. It was argued that 
if most of the galactic 'diffuse' background arises from such (unresolved) old SNRs then 
the absolute emissivity can be explained with an average magnetic field strength of (9(/xG), 
consistent with the RM measurements [28]. These large old SNRs may themselves generate 
the observed fluctuating component of the magnetic field. Given that they pervade the whole 
Galaxy and can be found at distances between hundreds of pc to tens of kpc, we expect them 



to contribute at angular scales between a few degrees up to the full sky — see Figs. 1 and 2. 





Figure 1. Sketch of the Galactic disk with shells 
of old supernova remnants and their skymap, i.e. 
their projection onto the celestial (hemi) sphere, in 
order to indicate the origin of 'Loop' features in the 
radio sky. The (helio-centric) Galactic coordinate 
system is shown by the grey lines. 



Figure 2. Zoom into the inner few kiloparsecs of 
the Galactic disk around the position of the Sun. 
The nearest old SNR which is responsible for the 
largest 'Loop' feature on the radio sky is circled. 



In this paper, we attempt to model the radio sky on a wide range of scales, taking 
into account the large-scale contribution of CR electrons and of the GMF, the variation of 
the GMF due to the turbulence cascade in the ISM and the effect of the presence of large 
radiative shells of old SNRs. The observable we model is not the all-sky radio map itself 
but its APS i.e. the distribution of Ci with /. Modelhng the APS instead of sky maps has 
several advantages; firstly, the cascading of energy from the outer scale L to the dissipation 
scale as well as the spatial distribution of old SNRs is a stochastic process of which the 
observed Galaxy represents only one particular realisation. What can be predicted are the 
ensemble averages which may to a certain degree differ from the observed Galaxy. The APS 
is a representation of this stochastic process that reduces the variation between different 
representations by efficiently averaging over many lines of sight. Secondly, the APS presents 
a representation of the information sorted by angular size and also encodes the symmetry 
properties of the sky map. Finally, for CMB studies (and also a particular class of indirect 
searches for DM annihilation signals), the observable is also the APS and it is therefore 
favourable to perform the foreground subtraction not on sky maps but directly on the APS. 

The remainder of this paper is organised as follows. In § 2 we discuss the APS of the 
'Haslam' 408 MHz all-sky survey [29]. In § 3 we present our model for the synchrotron 
background, discussing in detail each component (large-scale smooth emission, small-scale 
turbulence, intermediate scale old SNR shells) in turn. We demonstrate a good match of 
our model to the observed APS. In subsequent papers we will demonstrate a match to the 
observed frequency spectrum for different regions of the sky, as well as to the RM sky maps. 



2 The observed angular power spectrum 

We begin by computing the APS of the 408 MHz ah-sky survey [29]. Throughout this study, 
we used a verison of the map that has been corrected for scanning artefacts and with strong 
point sources removed.^ It is possible that unresolved point sources still contribute (in par- 
ticular at high multipoles), however we will argue later that this is most likely negligible. 
Furthermore, as this is the only full-sky point-source subtracted map that is publicly avail- 
able, we restrict ourselves, for the time being, to this one frequency. Important additional 
information is encoded in the frequency-dependence of the sky maps (or equivalently the 
APS), because of the different (and spatially varying) frequency dependence of the various 
physical components contributing to the diffuse emission, see § 3. This is particularly true 
for synchrotron emission by relativistic electrons, given that the spatial distribution of the 
radiating CR electrons varies with energy due to their energy-dependent diffusion and energy 
loss rates. Extrapolating the APS (or sky maps, for that matter) between widely separated 
frequencies is therefore non-trivial. (See [30] for an illustration of possible systematic ef- 
fects.) However, if successful in developing a convincing physical model for the APS, we 
can be of course more confident in predicting the APS at other frequencies. It would there- 
fore be desirable to test and (if necessary) calibrate our model by comparing to the APS at 
different frequencies. Here we adopt the complementary approach of fixing all the parame- 
ters using other observations, which allows for testing and dissecting the underlying (model) 
assumptions. 

In Fig. 3, we show the APS of the 408 MHz all-sky survey calculated using anaf ast from 
the HEALPix-suite [31]. For cosmological studies, it is customary to plot the Ci multiplied by 
/(/ + l)/(27r) in order to facilitate comparison with a Harrison-Zeldovich primordial spectrum 
for which C^ oc l//(/ + 1). Doing so is still helpful in the present context since multiplying C/ 
by /(/ + 1) ^ 1(21 + 1) directly displays the power contained in logarithmic intervals of /. 

It is seen that most of the power in the 408MHzmap is present at multipoles / < 10, 
corresponding to angular scales ^ 20°. The most prominent feature in the APS, however, is 
the disparity between even and odd modes — a consequence of the (approximate) symmetry 
properties of the synchrotron map, in particular the symmetry between the northern and 
southern (Galactic) sky and the invariance with respect to shifts in longitude: If a sky map 
J{0,(j)) is perfectly symmetric with respect to the Galactic plane, J{0,(j)) = J(7r — 0,0), 
then since the spherical harmonics satisfy Yij^^tt — 9,(f)) = (— l)^+^Y/^(0, 0), it follows that 
(^im — {~^y^^(^lm- If? in addition, we assume that the synchrotron flux is independent 
of longitude 0, then all a/^ vanish except when tti = 0. This results in the ai^ for even 
/ (and 771 = 0) being the only non- vanishing components, and therefore also the C/s with 
even /. However, as the above assumptions are satisfied only approximately (in particular, 
invariance with respect to shifts in longitude), the symmetry is partially broken and the C^s 
with odd / become populated. Thus the ratio of odd to even Qs is a measure of how much 
the Galactic plane symmetry is broken, and by how much the flux varies with longitude. 

It is also apparent from Fig. 3 that this odd-even disparity does not extend to arbi- 
trarily high multipoles. In fact, between / = 30 and 40, the saw-tooth structure largely 
disappears, pointing at a different component dominating on smaller scales. Moreover the 
APS at large (/ > 50) multipoles becomes rather smooth; this might seem surprising given 
how apparently stochastic structure is present in the sky maps on small scales. However, as 
mentioned already, the APS is formed at large / by averaging over many different m-modes 
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and correspondingly many different lines-of-sight (LoS), which leads to an efficient averaging 
of the stochastic information (i.e. the distribution of turbulent eddies or shells of old SNRs) 
and carves out the physically meaningful information that can actually be predicted, viz. 
their statistical expectation values. 
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Figure 3. The APS of the 408 MHz all-sky survey, both before (dotted line) and after (full line) 
correction for scanning artefacts and subtraction of point sources. 



3 Modelling of the angular power spectrum 

3.1 Definitions 

The power emitted at frequency v per unit solid angle and frequency by an electron of energy 
£^ in a magnetic field B is [32]: 



£(l/,^,B) 






(3.1) 



where B^ is the component of B perpendicular to the LoS, e, vn^ and 7 are, respectively, 
the electron charge, mass and relativistic gamma factor and i^5/3 is the Bessel function of 
degree 5/3. Given a number ne(r, E) of electrons per unit energy and volume at position r 
and with energy £", the total synchrotron power e{y^ r) is then 



z[y, t) = JdE n,{E, T)e{v, E, B(r)). 



(3.2) 



For a power law electron spectrum ne{E^r) oc £" ^, ^(i^, r) oc z/ ^^ ^^/'^B^^ ^^ (r). For the 
locally measured 7 — 3, this gives: ^(i^, r) oc Ty~^B'j_(r). 



The flux J{-u]9^(j)) observed in a particular direction (0,0) on the sky is obtained by 
integrating the synchrotron emissivity ^(z/, r) over the corresponding LoS, r(5, 0, (/>), 

J{i^;9,(t>)^ / dsi{,^,r{s,9,cf>))^ / ds dEne{E,r{s, 9, (f>))e{i^,E,B{r{s,e, ((>))) . 
Jo Jo J 

(3.3) 

In the following, we model the APS of the Galactic, diffuse radio background, taking 

into account dynamics over a wide physical range: the emissivity of the ISM with variations 

on large (see Fig 4a) and small (see Fig 4b) scales as well as the shells of old supernova 

remnants (see Fig 4 c). We also take into account free- free emission as well as the presence 

of unsubtracted point sources although this is not the main thrust of this work. 




(a) (b) (c) 

Figure 4. Different setups for calculations: (a) Smooth, large-scale emissivity with an off-center 
observer as considered in the GALPROP code (§ 3.2); (b) Turbulent halo dressed with the large-scale 
emissivity of the GALPROP code as in our MC calculation (§ 3.3); (c) As in (b) but with an additional 
contribution arising from a large number of old SNR shells (§ 3.6). 



3.2 Large-scale variation of synchrotron emissivity 

We start by modelling the GMF on scales much larger than L ^ 0.1 kpc, the assumed 
outer scale of turbulence. For this we use the GALPROP code which numerically solves the 
transport equation for Galactic cosmic rays (GCRs), taking into account the distribution of 
sources (as determined from e.g. SNR or pulsar surveys), the spatial diffusion of CRs in 
the turbulent GMF, their convection by Galactic winds, various energy loss processes and 
'reacceleration' due to diffusion in momentum space. We have adopted the parameter values 
of the "pure diffusion model" z04LMPD [33], which is known to reproduce both the locally 
measured electron spectrum as well as the synchrotron spectrum averaged over large parts 
of the sky out of the Galactic plane. These parameter values can also reproduce locally 
measured GCR fluxes and secondary-to-primary ratios. We note that the choice 6 = 0.5 in 
this model (z04LMPD) is somewhat at odds with the assumed Kolmogorv spectrum for the 
ISM turbulence (which would nominally predict 6 = 1/3). However, as a change of 6 can 
easily be absorbed into the source spectral indices without changing the primary spectra too 
much, we chose to retain the original values of z04LMPD. Note that as the assumed source 
distribution varies over distances > L, so will the electron density. Rather than describe the 
turbulent GMF fully, only its RMS value, averaged over ensembles of turbulent GMFs, is 
considered. For the normalisation and spatial dependence we have adopted the model that 



is stated [34] to give a good fit to global angular profiles of the 408 MHzall-sky survey, viz. 

B,^,ir,z) = Boe-^/''-\'\/^, (3.4) 

with Bq = 7.5 /xG, p = 30 kpc and ^ = 4 kpc. Note that we have later changed p to 100 kpc, 
i.e. essentially no variation with radius, to better reproduce the symmetry properties of 
the APS. We chose to ignore the ordered component altogether (for the time being) as its 
magnitude is much smaller than the RMS of the turbulent component. 

It has been recognised for some time [32] that the spectrum of the diffuse radio back- 
ground softens between at few hundreds of MHz and a few GHz which probably reflects a 
break in the power-law describing the underlying electron spectrum. We follow [34, 35] in 
adopting a spectral break in the (source) electron spectrum at T^i = 4 GV (which has been 
also found to be necessary for fitting the spectra of protons and nuclei [8, 36]). This break 
has also been attributed to a change in the energy dependence of the diffusion co-efficient 
[37, 38] but as the effect on the propagated spectrum is similar, we choose to remain agnostic 
as to its origin. The second break at 7^2 = 50 GV is probably not due to either cause but 
possibly the effect of an additional, harder contribution to the electron flux. (If the sources 
also accelerate positrons with a similar spectrum, this would explain the observed rise in the 
positron fraction [39].) The parameters adopted are shown in Table 1. 

The synchrotron sky map is readily obtained from the GALPROP code in the HEALPix 
scheme and its APS is again calculated using the HEALPix suite. 



Description 


Parameter 


Value 


Spatial diffusion coefficient at 4 GV 

Power lax index of the spatial diffusion coefficient 

Half-height of CR propagation volume 


^0 

5 

Zraax 


3.4 X lO^Scm^s"^ 

0.5 

4 kpc 


Source distribution 


pulsar-like 
for r > 10k 


40], but held const, 
pc 


for 7^ < 7^l 
Source spectral indices for TZi <TZ <TZ2 

for 7^2 < 7^ 

Break rigidities 


r2 

7^2 


1.6 
2.5 
2.2 
4GV 
50 GV 


Magnetic field strength at the Galactic center 

Scale height in radial direction 

Scale height in direction perpendicular to Galactic plane 


Bo 

p 


7.5 mG 

30 kpc, 100 kpc 

4 kpc 



Table 1. Parameters of the propagation and magnetic field model [33, 34] adopted for the computation 
of the large-scale component. 



3.3 Small-scale variation of synchrotron emissivity 

Clearly the assumption that the GMF has structure only on scales larger than L ^ 100 pc is 
unrealistic. Firstly the magnitude of the GMF as determined from Faraday rotation (which 
depends on the large-scale component, i.e. the GMF averaged over LoS > L) is too small 
[24] to reproduce the overall synchrotron flux (which depends on the mean square average 
value along the LoS). Hence, the variation on scales smaller than L cannot be neglected. 
Secondly, a turbulent GMF field is also inferred from the fact that the observed anisotropy in 



the arrival directions of GCRs is very small: around 10~^ up to 100 TeV and around 10~^ at 
lOPeV [41]. It is therefore believed that the dominant process of CR transport in the Galaxy 
is diffusion due to resonant wave-particle interactions — this requires a turbulent magnetic 
field on the scale of the particle gyroradius, e.g. ^ 10 pc for E = lOPeV in a /xG field. 

Under some simplifying assumptions the APS of the synchrotron emission from a turbu- 
lent magnetic field can be computed analytically [11]. Let the magnetic field be a Gaussian 
random field with an ensemble average that is only a function of the distance s from the 
observer: (B^) oc w{s) — we call this the statistically isotropic setup (see also Fig. 11). Let 
us approximate the electron spectrum as a power law n{E) oc E~^ which is close to the 
locally observed electron spectrum. The emissivity is then proportional to the square of the 
magnetic field perpendicular to the LoS which on average — assuming statistical isotropy 
— is 2/3 of the square of the total magnetic field. It turns out that the 3D power spectrum 
of the random field (i.e. the Fourier transform of the two-point correlation function) is then 
reflected in the APS of the resulting synchrotron sky map [11]. In fact, if we assume a Kol- 
mogorov power law with a cut-off at small wavenumber ko (corresponding to an outer scale 
of turbulence: L ^ l/feg), 

^ = J-2fc-^V3e-/cg//c^ . (3.5) 

If we further adopt w{s) = e~^ / ^ , then for / ^ 1, the APS is a broken power law in /: 

Ci = 2.96 X 10^ ^^^7;^ -^^ , (3.6) 



l + 5.0(^) J 



with power law indices —1 ( — 11/3) for / < /cr (^ ^ ^cr)- The break occurs around a critical 
multipole that depends on the path length R and the outer scale of turbulence L as /cr ^ 
{27TR/L) [11-13]. 

It has been argued that the observed APS follows a power law at large / with spectral 
index —11/3 due to this small-scale variation of synchrotron emissivity [11-13, 15]. However, 
for this to hold the normalisation of the observed high-/ component should also be reproduced. 
In fact, in the statistically isotropic setup, there is a close relation between the monopole Co 
and the power of the higher (/ > 0) multipoles. It is shown in Appendix A however that the 
power in the monopole is always many orders of magnitude larger than the power contributed 
by the rest of the APS. Therefore, fixing the normalisation Tq in eq. A. 4 to reproduce the 
APS at large / from eq. 3.6 will imply too large a monopole. We conclude that the magnetic 
turbulence can only contribute very little to the APS at large /. 

Of course the real question is how the picture changes when we consider a more realistic 
setup. We now consider an off-center observer in a medium with small-scale turbulence, that 
is dressed with the overall large-scale emissivity of the GALPROP run from § 3.2 (cf. Fig. 4b). 
One would then expect the APS to be a mixture of the large-scale contribution and the 
above small-scale APS. In fact the APS should be more-or-less the sum of the large-scale and 
the small-scale APS as these are disjoint scales. For the small-scale component, one would 
need to average the APS over a distribution of scale-lengths R since the scale length in the 
off-center setup depends on the direction [12]. Since the large-scale emissivity is no longer 
isotropic, the monopole computed in (A. 2) gets spread out into the lowest multipoles and we 
would therefore expect the large-scale power in these lowest multipoles to dominate over the 
small-scale power at the higher multipoles. 



We have computed the sky map and its APS numericahy by a Monte Carlo calculation. 
To that end we have prepared a turbulent magnetic field by superposing hundreds of plane 
waves evenly spaced logarithmically in k with amplitudes chosen so as to reproduce a power 
law power spectrum oc k~^^/^ for k > 2ti / L [42] having unit variance. In particular, we adopt 
a coherence length L ^ 400 pc. This statistically isotropic Gaussian random field b(r) is then 
considered as the deviation of the turbulent magnetic field B(r) from its ensemble average 
5rms(^5^) (eq. 3.4), i.e. B(r) = Srms(^5 ^)b(r). Given the (smooth) emissivity (cf. § 3.3) 
calculated by GALPROP, s{v]B^^Q{r^z))^ we can now compute the emissivity from the fully 
turbulent field by using the scaling of the emissivity and critical frequency with the magnetic 
field B (cf. eq. 3.1), 

5(^;/3S)=/35(^//3;S), (3.7) 

where /3 = |b^(r)|. We generate a sky map by integrating this emissivity along LoS's over 
the whole sky and calculate its APS as above. 

3.4 Free-free emission 

Free- free emission is the bremsstrahlung of thermal electrons in the ISM. For an electron 
temperature Tg (in units of degrees Kelvin), the optical depth is given by [43] 

Tc = 3.014 X lO-^T^^/^z/"^ (In [4.955 x lO-^z/"^] + 1.51nTe) / drnen^ (3.8) 

JLoS 

and the brightness temperature is oc TcTq. A reasonable approximation is 

Tc oc T-^-^^v-^-^ I dr nefii , (3.9) 

JloS 

and the brightness temperature then factorises into a frequency part (oc z/~^-^) and a spa- 
tial/angular part. 

Consequently, unlike synchrotron emission the spatial distribution of thermal brems- 
strahlung emissivity is not expected to change much with frequency (unless the ISM electron 
temperature varies significantly between locations). We can therefore attempt to trace the 
free-free emission with a template. Here, we employ the estimate for the free-free emission 
at 23 GHz provided by the WMAP team using a maximum entropy method (MEM) [44] and 
their adopted frequency spectral index of —2.15. 

We note that the MEM free-free template is not an unbiased tracer of thermal brems- 
strahlung. First, it was derived by separating the observed 23 GHz sky map into synchrotron, 
free-free and thermal dust emission and therefore any error in the synchrotron extrapolation 
(which we try to improve upon here) translates into errors on the MEM free-free template. 
Moreover, spinning dust emission was ignored all together such that the MEM templates po- 
tentially compensate for this. The template therefore likely overestimates the real free-free 
emission. Alternatives to the MEM template, e.g. the WMAP MCMC model, have been 
shown however to be even less reliable [21] and over- (under-)estimate free-free emission when 
excluding (including) a spinning dust component. The MEM free- free template on the other 
hand has a rough correlation with maps of radio recombination lines which is why we adopt 
it here. 
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3.5 Unsubtracted point sources 

It is possible that the point source subtractions performed on the 408 MHz ah-sky survey have 
not been able to subtract all the sources present in the sky map. Relatively faint sources, 
and even relatively strong sources in the Galactic disk, might have been missed. Here, we 
model their APS by shot noise, as is expected for an unclustered population of point-like 
objects (see [45]), and fix the normalisation such that the observed APS is matched at the 
highest / accessible, i.e. / ^ 200. We find that Q ^ 2 x 10~^ K^ allows for a good fit at large 
/. By comparing the APS of the unsubtracted with that of the (partly) subtracted 408 MHz 
all-sky survey, we have confirmed that this seems to indeed describe the APS of the point 
sources rather well. 

3.6 Shells of old SNRs 

As we will see in the next section, when the APS is modelled by adding contributions from 
both large- and small-scale variations of synchrotron emissivity, as well as free-free emission 
and unsubtracted point sources, there is still a large deficit in power with respect to the 
observed APS on intermediate scales (10 < / < 100) that cannot be explained either by an 
increase in free- free emission, nor the presence of unsubtracted point sources. However, a 
class of extended sources that would contribute to the APS on such scales is visible already in 
sky maps from radio to microwaves (and even in gamma-rays): the so-called 'radio loops' [46]. 
These ring-like structures, that span up to 60° in angular radius, are believed to be the large 
(up to ^ 200 pc in radius) shells of old SNRs (see [47] for a recent study of Loop I, also called 
the North Polar Spur). 

Towards the end of its lifetime, a SNR is in its radiative phase with a low expansion 
velocity of (9(100) kms~^ such that particle acceleration is no longer expected to take place. 
However, as the expansion is now in the pressure-driven 'snow plough' phase [48], the shock 
compression factor 77 can get very large leading to a strong enhancement of the ISM magnetic 
field upon crossing the shock. (In the following we adopt 77 = 6 which is rather conservative.) 
This has two important consequences [26, 28]: First, by betatron acceleration the energy 
of electrons increases by a factor ^ y^r/ upon crossing the shock; as the electron spectrum 
is a steeply falling power law, this increases the emissivity for a fixed frequency. Second, 
due to the compression of the magnetic field by the same factor 77, the emissivity at fixed 
frequency is now sourced by electrons which are a factor 77 lower in energy — these are more 
numerous, thus further increasing the emissivity. Hence, a compressed shell easily outshines 
the emissivity of the surrounding ISM. 

Given at set of shells at positions {r^}, the total synchrotron flux Jsheiis(^;^5 ^) is the 
sum over all individual shells i, and so are the aim'- 

Jshellsil^, i,b) = J2 Js^eiuil^^ ^' b) ^ «'-"' = E «?-''' • (3-10) 

i i 

We approximate the compressed medium of an old SNR as a spherical shell of inner 
(outer) radius Ri {R2) with a constant compression factor 77. Then the flux from a single 
shell i centred at at r^ factorises into a frequency- dependent specific emissivity ei{v) and an 
angular LoS integral gi{l^b): 

Jsheiui^J.b)^ J dEn^,,^a{E;v,)e{u,i^B{v,),E) J dsp{\T{s,9,^^ (3.11) 
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where p{r) = 1 for i?i < r < i?2 and otherwise. Here, we have approximated both the 
betatron- accelerated electron spectrum nE;a^c&d(^) and magnetic field B(r) in the shell by 
their values at its centre n, i.e. nacc'd (E^Vi) and B(r^). 

3.6.1 Specific emissivity of a shell 

Consider an electron with pitch angle Oq^ i.e. momentum p± = psin^o {p\\ — pcos^o) 
perpendicular (parallel) to the magnetic field. Due to efficient compression of the GMF, the 
magnetic field in the shell will be mostly aligned with the shock, that is tangential to the 
shell. If the shock width is larger than the gyroradius of the (GeV) electron, the time-scale 
over which the magnetic field B changes is larger than a gyration period, and therefore the 
first adiabatic invariant p\/ B is conserved. Upon crossing the shock with compression factor 
77, the electron therefore gains a factor ^Jr] in p^, while p\\ remains constant. An initially 
isotropic distribution n{p) therefore becomes anisotropic: 

n{p) ^ nacc'db, e)^n[ P | . (3.12) 

where the new pitch angle 6> is a function of the initial pitch angle Oq (see [26] for details). 
We assume that this distribution gets quickly isotropised by pitch-angle scattering through 
turbulence behind the shock. The distribution then needs to be averaged over the new pitch 
angle 6, however, an order of magnitude approximation is given by: 

^acc'd(p, 0) ^ — n vyY ) ' ^^'^^^ 

assuming that only the two perpendicular components of the momentum change. 

The emissivity ^acc'd ^tt frequency u of the compressed shell can now be simply expressed 
in terms of the nominal emissivity i of the (unaccelerated) electrons in the (uncompressed) 
magnetic field, but at a different frequency: 

4cc'd(^, VB) = JdE^n (eJ^J e{u, r]B, E) (3.14) 

= ^J^JdE'n{E') e L.r^B, \j\r^E'\ (3.15) 



^^^ismIt^,^) , (3.17) 



where we have used the scaling of the specific synchrotron emissivity (see eq. 3.1). 

3.6.2 Line-of-sight integral over a shell 

The projection of the shell onto the celestial sphere is a limb-brightened disk and for a given 
distance d^, the flux depends on the angular distance 9' — arccosz'' from the projected centre 
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of the sphere only (see Figs. 5 and 6): 



9i{z') = 2 



Rl - df (1 



yl2\ 



Rl 



dm 



y'2\ 



Rj-d^il-z''^) for 



djil-z'^)<Ri, 



for Ri< ^dj{l-z'^)<R2, 
otherwise . 

(3.18) 





radio profile of Loop I (Salter 1983) 

shell of const, emissivity (i?i = 0.7 d , R2 = O.S d ) 



20 30 40 50 

angular distance from centre of shell [°] 



Figure 6. Angular profile of a shell of an old 
SNR. The data points show the (frequency aver- 
aged) angular radio profile as measured in Loop 
I [49]. The solid line is the profile of a shell 
of constant volume emissivity with inner radius 
Rl = 0.7 d and outer radius R2 = 0.8 d. 



Figure 5. Coordinate systems used for the 
computation of the APS of the shells. The 
black, unprimed coordinate system O is helio- 
centric with the Galactic centre at (x, 7/, z) = 
(8.5,0,0)kpc, and the Galactic north pole in 
positive z-direction. The blue, primed coordi- 
nate system O' is centred on the direction of 
the centre of a shell and can be obtained by a 
rotation of ii around the z axis and a rotation 
of 7r/2 — bi around the 7/'-axis. 



We start by considering the a^^ of an individual shell z in a coordinate system O^ with 
z^'-axis in the direction of the centre of the shell (see Fig. 5), 



^Si{u)J dz' J 



2- ^J2l + 1 lil-m)\^^^,^j^.. 



d(P' 



47r 



(l + m] 



-Pr{z')e^'^^ g{z') 



^ii'')^J'-^I-idz'Piiz')g{z') for m = 0, 

otherwise . 



(3.19) 
(3.20) 
(3.21) 



The spherical harmonics coefficients a) in galactic coordinates O can be computed 
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with the Wigner D-matrices Df"^{-£i, 7r/2 - h, 0), 



air 



i{v) J2 A'"''"(-^i,|-^i,0)aL 



m'=-l 



•\0m 



-£i 



TT 



z,7r -h,0 



21 + 1 



An 



= ei{v)D\ 

= ei{v)yL (I - hu -^^) j_^ dz'Pi{z')g{z') . 



l'^dz'Pi{z')g{z') 



(3.22) 

(3.23) 
(3.24) 



The integral in the last line can be written: 



J dz'Piiz')g{z') = dJii^l - {R2/dY) - dJii^l - {Ri/df) , (3.25) 



where 



Jl{a)^ f dz'Pi{z')Vz^-a^ 

J a 



(3.26) 



can be solved analytically in terms of Gauss hypergeometric functions 2F1 (cf. App. B). 

3.6.3 The Galactic population of radio loops 

Only four of the radio loops are visible in the sky and can be modelled directly. We fix their 
radius to 200 pc which is consistent both with calculations of the evolution of SNRs [48, 50] 
and estimates for Loop I [47, 49]. We have adopted the positions and angular sizes as found 
originally [46] and computed the distances, see Table 2. 



Object 


I (centre) 


b (centre) 


diameter 


distance 


Loop I 


329° 


17.5° 


116° 


240 pc 


Loop II 


100° 


-32.5° 


91° 


280 pc 


Loop III 


124° 


15.5° 


65° 


370 pc 


Loop IV 


315° 


48.5° 


39.5° 


590 pc 



Table 2. Angular positions, diameters and resulting physical distances of the local radio spurs for an 
assumed physical radius of 200 pc [46] . 

However, the observable loops constitute most like only a small fraction of the total 
number of old SNRs in the Galaxy. Adopting a SN rate of 1 — 3 x 10~^ yr~^ and estimating 
the duration of the radiative phase to be between 10^ — 10^ yr, we can expect up to thousands 
of those old SNR shells to contribute to the synchrotron emission in radio and microwaves. 
Most of them will however not be directly observable in the sky since being further away than 
the 4 local loops they will blend into the diffuse emission from the Milky Way. We therefore 
model them by Monte Carlo by randomly choosing (9(1000) positions distributed according 
to the pulsar-like source distribution adopted by GALPROP (see Sec. 3.2 and Table 1). As 
we will add the contribution from this random set to that from the known local loops, we 
have excised from the random set all those shells that are closer than 500 pc in order not to 
double-count the local shells. 

We fix the radius of all shells to 200 pc and calculate their spherical harmonics expansion 
coefficients aj^ from eqs. 3.17 and 3.25. Finally, we sum up the spherical harmonics expansion 
coefficients to a^^^^^ and compute their angular power spectrum Ci^^^^^. 
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4 Results 

4.1 APS from conventional model 
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Figure 7. Top panel: The APS of the synchrotron sky computed by taking into account the 
large-scale and small-scale variations of emissivity (orange line) , free-free emission (purple line) and 
unsubtracted sources (grey line). The sum of these components (blue line) is compared with the APS 
from the 408 MHz all-sky survey (black line). Bottom panel: The relative residual of the modelled 
APS with respect to the observed one. 

In Fig. 7 we show our modelling of the angular power spectrum resulting from both 
large-scale and small-scale variations of synchrotron emissivity as well as free-free emission 
and unsubtracted point sources and compare it to the observed APS of the 408 MHz all-sky 
survey. For all modelled components we have taken into account the exponential cut-off due 
to the beam window function. 

We start the discussion by computing the APS from the variation of the synchrotron 
emissivity on large scales only (orange line). We note that most of the power is indeed 
contained on large scales (/ < 10) and the APS falls off quite rapidly for higher /, as is 
indeed expected given that all the ingredients (source distribution, transport parameters and 
magnetic fields) vary only on kpc scales. The APS computed with the parameters of [33, 34] 
does show the expected even-odd multipole structure but the difference between the power 
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in even and odd multipoles is too small. This signals that the distribution in longitude is not 
extended enough; if it was more extended this would bring us closer to the symmetric case 
and would suppress power in the even multipoles. In Fig. 7 and henceforth we have therefore 
adopted p — 100 kpc. 

Next we compute the contribution of the ISM synchrotron emissivity, with (red line) 
and without (orange line) the small-scale turbulence. We note that while the turbulent APS 
indeed shows a power law behaviour C^ oc / for / > 50, the additional power due to turbulence 
is relatively small. This corroborates our arguments from § 3.3 about the relative power in 
the monopole and in higher multipoles and demonstrates that the large-scale behaviour of 
the measured APS cannot be due to ISM turbulence. 

Now adding the contribution from free-free emission (yellow line) and the Ci due to 
unsubtracted sources (grey line, almost negligible below / ^ 100), we arrive at the total 
APS (blue line). Comparing this to the measured APS, it is evident that ISM synchrotron 
and free-free emission alone can give a reasonable fit only for the dipole and quadrupole 
and around / ^ 100. At all other /, and in particular at intermediate multipoles, the model 
underproduces the power observed in the 408 MHz all-sky survey. 

4.2 APS adding SNR shells 

Figure 8 shows the modelled APS as before, but now including the contribution from the 
four local plus the ~ 1000 shells from the Galactic random distribution, for 10 different 
realisations of this ensemble. The APS for the shells alone is also shown and it is clear that 
the shells mostly contribute to the total APS on large and intermediate scale, i.e. / < 50 
which is where we previously found the deficit. The agreement with the measured APS is 
now much improved and the residuals are at most 20%. We note that the shells do not 
dominate the powerfor any / so the slightly larger residuals for / > 100 must be attributed 
to the details of the free-free MEM map. 

It is also note-worthy that the differences between different realisations of the distri- 
bution of Galactic SNR shells is relatively small. This illustrates our point about the APS 
being particularly appropriate for modelling of Galactic foregrounds as the influence on the 
APS of unknown parameters (e.g. the position of and distances to the shells) is relatively 
mild. 

In Fig. 9 we show the APS again including all contributions and for the best fit to the 
observed APS. Our model reproduces the observed APS with residuals of order 10% up to 
/ ^ 200, i.e. down to a degree. 

5 Discussion 

The reader might wonder why we attempt to develop a physical model of the unpolarised 
synchrotron emission in order to e.g. do a reliable foreground subtraction for studies of 
the CMB temperature anisotropics, given that it is apparently possible to study the CMB 
temperature anisotropics rather successfully without requiring a fully physical model of the 
foregrounds [51]. This is justified by saying that most foregrounds (in particular on small 
scales, i.e. at high /) are contained in the Galactic disk and can therefore be suppressed by 
simply masking this part of the sky. However apart from complicating the analysis (spherical 
harmonics on the cut sky are not orthonormal, thus introducing correlations between the aim 
of different /), it is hard to estimate by how much the masked sky maps are still contaminated 
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Figure 8. Top panel: The APS of the synchrotron sky computed by summing the large-scale and 
small-scale variations of emissivity, free- free emission, unsubtracted sources as well as the contribution 
from different configurations of shells of old SNRs (also shown separately). Bottom panel: The 
relative residual of the modelled APS with respect to the observed one. 



by high-latitude Galactic foregrounds. It is quite possible that the low multipole anomalies 
of the CMB are at least partly due to such residual contamination. 

Second, the modelling of the unpolarised synchrotron emission is of course only a first 
step, with the next frontier being the observation of B-mode polarisation which will necessi- 
tate a sound understanding (if not a reliable subtraction) of the polarised emission. In fact, 
the expected signal for realistic r is only about 12% of the foreground on degree scales over 
75% of the sky, hence to reliably extract the signal would require an understanding of the 
polarised foreground down to the few percent level [3]. While our fit is admittedly accurate 
only to 10 % or so, we feel that — compared to previous studies of the synchrotron emission in 
sky maps (as opposed to our study of the APS) — this already a considerable improvement. 
We emphasise that this has not been achieved by means of a detailed parameter study of 
existing synchrotron models but by an alternative representation — the APS — of the infor- 
mation contained in the radio sky map, and comparison with the simplest possible models, 
i.e. the standard GALPROP setup and simple physics describing the shells of old SNRs. We 
believe that this shows the power of the APS approach for modelling diffuse foregrounds — 
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Figure 9. Top panel: The APS of the synchrotron sky computed by summing the large-scale and 
small-scale variations of emissivity, free- free emission, unsubtracted sources as well as the contribution 
from the best fit configuration of shells of old SNRs (also shown separately). Bottom panel: The 
relative residual of the modelled APS with respect to the observed one. 



applications are certainly not constrained to radio and microwave frequencies alone. 

As to the changes necessary for modelling the polarised emission, one needs of course 
to model the ordered component of the GMF, but in principle the other ingredients remain 
the same. In particular if we follow the picture of old SNR shells as compressions of the 
ISM together with the GMF, the magnetic fields should be mostly tangential which should 
simplify the modelling of the polarised emission from the shells. This simple picture seems to 
be borne out by the polarised emission observed from the four local loops in the microwave 
range, with polarisation vectors mostly radially orientated (see Fig. 10 of [52]). Furthermore, 
given the lack of free-free emission in polarisation, we are confident that our approach will 
allow a similar (if not better) level of precision in polarisation. 

Although we have stressed the validity of modelling of the APS alone, for completeness 
we also show in Fig. 10 sky maps of the different components contributing to the diffuse 
emission at 408 MHz. We note that with the inclusion of turbulence up to an outer scale 
of 400 pc, the diffuse emission of the ISM (top left panel) looks very similar, in particular 
the level of "patchiness" is similar. However it appears that the emission in the Galactic 
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total map at 408 MHz 



Haslam 409 MHz 




Figure 10. Top left: ISM synchrotron emission on large (modelled by GALPROP) and small 
scales. Top right: Free-free emission, extrapolated from the WMAP MEM free- free map. Middle 
left: Synchrotron emission from the Galactic population of shells of old SNRs. Middle right: 
Synchrotron emission from the four local shells of old SNRs, viz. Loop I - IV. Bottom left: The sum 
of synchrotron emission on large and small scales, free-free emission. Galactic population of shells and 
local shells. Bottom right: The 408 MHz all-sky survey. 



plane is not as peaked as is observed (lower right panel). The contribution from the Galactic 
population of shells (middle left panel) on the other hand does contribute in the plane, 
given that the majority of the shells is very far away, and brings the latitudinal profile into 
agreement with the observed one. Note that there are still a couple of shells at distances 
< 1 kpc that lead to visible loops in the sky map. The most obvious difference to the sky 
map for the local, i.e. those with distances < 500 pc, shells (middle left panel) apart from the 
larger angular size, is their offset from the Galactic plane. (We have constrained the centres 
of the non-local shells in our MC simulation to lie in the Galactic plane but loosening this 
constraint does not change the APS much.) Finally, the sum of all components (bottom left 
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panel) shows a remarkable (qualitative) agreement with the observed sky map (bottom right 
panel). We note that the observed radio loops do not possess the perfect ring-like structures 
as assumed in our simplified approach in § 3.6; the emissivity probably varies due to the non- 
uniform evolution of the SNR, the overall structure of the ordered GMF and the interaction 
and overlap of shells. However, in the APS this does not seem to lead to too big a deviation, 
as it is mostly sensitive to the overall angular structure. 

6 Summary 

We have presented a new model for the diffuse. Galactic synchrotron radio background that 
takes into account variations of the emissivity from the scale of the Galaxy as a whole, down 
to the scale of interstellar magnetic turbulence. We uncovered a lack of power on intermediate 
scales when only the variation of the ISM emissivity is taken into account (besides free-free 
emission and possibly unsubtracted point sources). We argued that this deficit is most likely 
cured by the inclusion of angular correlations due to the presence of thousands of old SNR 
shells (similar to the known four local 'radio loops'). Even in the simplest approach, i.e. 
when we fix the free parameters using other observables, the agreement with observations 
at 408 MHz is quite remarkable. This is mainly due to the chosen observable, i.e. the APS, 
being particularly appropriate in the sense that it averages out much of the uninteresting 
stochasticity and highlights the dependence on the underlying physically interesting model. 
We believe that the agreement can be improved even further by calibrating our model to 
point source subtracted all sky maps at other frequencies, for example at 1.42 GHz [53-55], 
which have so far not been publicly available. 
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A The monopole compared to the total power 

Given a (rescaled) synchrotron skymap, 

j{e, 4>) = j ds w{s)Bl{r{s, e, (/.)), (A.l) 

the monopole is proportional to the squared average of this sky map: 



The total power is, on the other hand, defined as the square of the map integrated over the 
whole sky, 

Ctot= / df^J*(0, 0)7(0, 0) = 5^(2/ + 1)C/, (A.3) 

i.e. the weighted sum of the APS. In the present case, i.e. for a power law turbulence 
spectrum 

^^'=J-2fc-ii/3e-'=eA% (A.4) 



d3fc 
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it can be computed via eq. 3.6. In Fig. 12 we show (Ctot — Co)/Co, i.e. the power in modes 
/ > relative to the power in the monopole, as a function of L/R. This is essentiahy the 
variance of a sky map J(0, 0) divided by its mean square: 



Ctot - Co ( J2) - {Jf 



Co 



{J? 



(A.5) 



where the angled brackets denote averages over the sky. It is apparent that the multipoles 
/ > only contribute to the total power of the sky map when L/R > 1. However, this would 
mean that the outer scale of turbulence is larger than the column depth - a contradiction in 
terms. We conclude that for realistic parameter combinations, e.g. L ^ 100 pc and R > 1 kpc 
the monopole always dominates the total power of the map. This is essentially a consequence 
of the central limit theorem: as the number of turbulent eddies that is averaged over by the 
LoS integration becomes large, the variance around the average gets small. 
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Figure 11. Observer at the centre 
of a turbulent sphere as in the an- 
alytical computation for small-scale 
varying turbulence. 



Figure 12. Power in modes / > relative to the power in 
the monopole, as a function of LjR. 



B Auxiliary function J/ (a) 

We start evaluating the integral J/ (a) by partial integration and by using recurrence relations 
of the Legendre polynomials. 



Ji{a) = \ dz'Pi{z')^/^^. 

J a 



^(P;+l(z)-P;_l(z))Vz2_„2 



Vz^ — a 



21 + 1 Ja ^/^^^^\2l + 3 



2-a2 

dz / 1 + 2 



iPi+iiz) - Pi.iiz)) 



Pl+2{Z) + 



I 
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2^ + 3 21-1 



Pl{z) 



1-1 
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(B.l) 

(B.2) 

(B.3) 

Pl-2{z)] . 



(B.4) 



21 



Integrals of this form compute as the real part of the same integrand but with the lower 
boundary at 0, 



/ 

J a 



.'v^^^'^^U't^^*'^- <^^^' 



^z^ — o? 



The integral in the brackets can now be solved by exploiting that d(arsinh z)/(\z — l/\/l + z^ 
and expanding arsinh z into a power series. The Legendre transform is evaluated separately 
for odd and even powers in z but in both cases, ones finds 



-/" 

<^ Jo 



^z _ / X a/tt ~ / 1 1 13 1 1 . /_ X 

Pliz) = i-3F2 1, o' o! 1 - o' o + o; 9 ' (B-6) 



.^x2 '^ ^ 2a'' ^ V '2'2' 2'2 2' a2 



where we have set a = ia and 3F2 is the regularised hypergeonietric function. After some 
algebra, this reduces to 

/ Pi{z) = ^^-^ ^Fi{- + -,- + -;- + l; . . (B.7) 



a Jo /(L)^^i 20F ^ ' V2 2' 2 2' 2 ' a2 

Taking into account that a is purely imaginary, we finally find the real part 

(B.8 
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